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. We report on the status of the calculation of deep-inelastic structure 

' functions at three loops in perturbative QCD. The method employed allows 

, to calculate the Mellin moments of structure functions analytically as a 

f^) ■ general function of N. As an illustration, we present the leading fermionic 

contributions to the non-singlet anomalous dimension of F2 at three loops 
and, as a new result, to the non-singlet coefficient function of F% at three 
loops. 
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1. Introduction 



Today, structure functions in inclusive deep-inelastic scattering are ex- 
tremely well measured quantities. As a consequence, they offer the possibil- 
ity for very precise determinations of the strong coupling a s and the parton 
distribution functions. The high statistical accuracy of the present and 
upcoming experimental measurements demands analyses to next-to-next- 
to leading order (NNLO) of perturbative QCD for the structure functions 
F 2 ,F 3 andF L . 

However, the complete NNLO corrections are not fully available yet. The 
two-loop coefficient functions of i^j-^b and Fl have been calculated [1, 2]. 
For the three-loop anomalous dimensions only a finite number of fixed Mellin 
moments [3] are presently available. In addition, some information about 
leading fermionic contributions [4, 5] and the small- x limit [6] exists. 

In the following, we briefly report on the status of the calculation of 
the coefficient functions and the anomalous dimensions to three loops in 
perturbative QCD. Furthermore, we present results for the leading fermionic 
contributions to the non-singlet structure function F^. 
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2. Method 

We employ the optical theorem and the operator product expansion 
(OPE) to calculate the deep-inelastic structure functions in Mellin space 
analytically [7, 8] as general functions of N. For the N-th Mellin moment 
of F2 we can write 

F 2 N (Q 2 ) = J dxx N - 2 F 2 (x,Q 2 ) = ]T (%,a s )Ai N (^) , (1) 

j=a,q,g V V I 

where C|j denote the coefficient functions and A J P N the spin averaged 
hadronic matrix elements of singlet operators O q , O g and non-singlet op- 
erators O a , a = 1, 2, ... , (n 2 — 1), of leading twist. Both, the coefficient 
functions and the renormalized operator matrix elements in eq.(l) satisfy 
renormalization group equations governed by the same anomalous dimen- 
sions 7jfc. The anomalous dimensions determine the scale evolution of deep- 
inelastic structure functions. 

The calculation of the coefficient functions C^j and anomalous dimen- 
sions 7^ at a given order in perturbation theory amounts to the deter- 
mination of the iV-th moment of all contributing Feynman diagrams with 
external partons of momentum p with p 2 = and photons of momentum q 
with q 2 = —Q 2 . To achieve this task, we apply the following strategy [2, 9]. 
We set up a hierarchy among all diagrams depending on the number of p- 
dependent propagators. We define basic building blocks (BBB) as diagrams 
in which the parton momentum p flows only through a single line in the 
diagram, while composite building blocks (CBB) denote all diagrams with 
more than one p-dependent propagator. 

Then, with the help of integration-by-parts [10] and scaling identities [2] 
we determine reduction schemes that map the CBB's of a given topology 
to the BBB's of the same topology or to the CBB's of a simpler topology. 
Subsequently, we use reduction identities that express the BBB's of a given 
topology in terms of simpler topologies. Working in Mellin space, the re- 
duction equations often involve explicitly the parameter N of the Mellin 
moment. Sometimes, one encounters difference equations in N for the iV-th 
moment F(N) of a diagram, 

a (N)F(N) + ai (N)F(N - 1) + (2) 
... + a n (N)F(N - n ) + G(N) = 0, 

where G(N) denotes the iV-th Mellin moment of simpler diagrams. First 
order difference equations can be solved at the cost of one sum over T- 
functions in dimensional regularization. We use D = 4 — 2e. The T- 
functions can be expanded in e and the sum can be solved to any order 
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in e in terms of harmonic sums [11, 12]. Higher order difference equations 
could be solved constructively. On the mathematical side, the approach to 
calculate Mellin moments of structure functions relies on particular mathe- 
matical concepts [13], such as harmonic sums [11, 12] and our ability to set 
up and solve the difference equations as nested sums in N. 



3. Leading fermionic contributions 

To illustrate the method, we discuss the leading fermionic contributions 
to the non-singlet structure function. At three loops, they are proportional 
to np with rif being the number of massless fermions. These contributions 
form a gauge- invariant subset, but do not involve yet any genuine three- loop 
topologies. Therefore, in the sense of the reduction strategy sketched above, 
the n^-terms are easier to calculate. 

The result for the n^-contribution to the non-singlet anomalous dimen- 
sion <7qq^' ns at three loops is known from the work of Gracey [4]. It is given 
by, 

(2)ns „ 2 A 7 32 1 88 1 8 1 32 1 
g (z),ns _ (j Z 1 1 

yqq F f \9 ^ 9 N+l 27(Af+l) 2 9(Af+l) 3 9 N 

with C F = (N 2 - 1)/(2N), which is 4/3 for QCD. 

As a new result, we give here the n^-contribution to the non-singlet 

coefficient function c^™ s at three loops for the flavour class, where both 
photons couple to the external quark. Strictly speaking, the three-loop 
coefficient functions contribute in a perturbative expansion only at next- 
to-next-to-next-to- leading order (NNNLO). However, the result illustrates 
nicely that our method will not only provide the anomalous dimensions, 
which are proportional to the single pole in e in dimensional regularization, 
but also the coefficient functions which are determined by the finite terms, 
since our approach is (at least in principle) not limitated to a given order 
in e. 

In terms of harmonic sums up to weight four, our result for the n 2 - 
contribution to ci 3 ^' ns reads 

(3), ns _ 



C 2,qq 



_ 2/ 9517 8 . 3 6 7 48 1 16 C3 4384 

CFn f —AM ~ Q C3 + 790 AfZl + 



486 9^ 729 N+l 27N+1 81 (N+l) 5 
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2360 1 184 1 16 5i(iV+l) 544 5i(iV+l) 

"I 77, / TV T . 1 \ Q /AT . -I \ A 



81 (N+l) 3 27 (N+l) 4 3 (N+l) 3 27 (N+l) 2 
32Si,i(iV+l) 32 5 2 (N+1) 2240 5i (N+l) 272 Si,i(JV+l) 
~9~ (N+l) 2 + ¥ (N+l) 2 + ~81 N+l + ~2T N+l 
16 5i,i,i(JV+l) 16 5i, 2 (N+l) 272 5 2 (N+1) 16 S 2 ,i{N+l) 
9 N+l 9 N+l 27 N+l 9 N+l 

^16 5 3 (N+1) 11170 1 16C 3 + 1204 1 992 1 184 1 



+ 



N+l 729 AT 27 N 81 N 2 81 N 3 27 N 4 

16 Si(N) 232 Si{N) 32 S ltl (N) 32 5 2 (N) 644 gi(N) 

~Y N 3 + "27" N 2 + ¥ N 2 T^V 2 "81 N~ 
104 5i,i(7V) _ 16 51,1,1 (JV) 16 5i, 2 (7V) 104 5 2 (N) 

27 JV 9 N 9 N 27 N 

l 16 5 2 ,i W 16 5 3 (N) 8714 ,32 940 

232 32 32 232 

+^rSi,i,i(iV) + -5i,i,i,i(JV) - ¥ 5i,i, 2 (iV) - — 5i, 2 (iV) 

-f Si,2,i(N) + f S lfl {N) - ^S 2 (N) - ^S 2 ,i(iV) - yS 2 ,i,i(iV) 

+ Y S2 ' 2(N) + lt S3iN) + 1 5s ' l(iV) " ^T S ^ N) ) ■ (4) 



This result agrees with the one of the fixed Mellin moment calculation [3] 
for N = 2, ... ,14. It may be used to improve approximations [14] to the 
full functional form of c^™, although the n 2 -terms are numerically not the 
most dominant contribution. 



4. Conclusions 

The present approach based on the OPE, to calculate the Mellin mo- 
ments of structure functions allows for the calculation of the complete co- 
efficient functions and the anomalous dimensions at three loops. 

As a next step, one can consider the subleading fermionic contributions 
at three loops proportional ton/. The determination of these terms already 
relies on major parts of the complete reduction scheme as it requires the 
calculation of several genuine three-loop topologies of the Benz type as well 
as the calculation of two-loop topologies with a self-energy insertion. The 
results for the n /-contribution to the non-singlet anomalous dimension will 
be presented elsewhere. 
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